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ABSTRACT 

Context. The properties of variable stars can give independent constraints on the star formation history of the host galaxy, by 
determining the age and metallicity of the parent population. 

Aims. We investigate the pulsation properties of 84 Anomalous Cepheids (ACs) detected by the OGLE-III survey in the Large 
Magellanic Cloud (LMC), in order to understand the formation mechanism and the characteristics of the parent population they came 
from. 

Methods. We used an updated theoretical pulsation scenario to derive the mass and the pulsation mode of each AC in the sample. 
We also used a Kolmogorov-Smimov test to analyze the spatial distribution of the ACs, in comparison with that of other groups of 
variable stars, and connect their properties with the star formation history of the LMC. 

Results. We find that the mean mass of ACs is 1.2 + 0.2Mo. We show that ACs do not follow the same spatial distribution of 
classical Cepheids. This and the difference in their period-luminosity relations further support the hypothesis that ACs are not the 
extension to low luminosity of classical Cepheids. The spatial distribution of ACs is also different from that of bona-fide tracers of 
the old population, such as RR Lyrae stars and population II Cepheids. We therefore suggest that the majority of ACs in the LMC 
are made of intermediate-age (l-6Gyr), metal-poor single stars. Finally, we investigate the relation between the frequency of ACs 
and the luminosity of the host galaxy, disclosing that purely old systems follow a very tight relation and that galaxies with strong 
intermediate-age and young star formation tend to have an excess of ACs, in agreement with their hosting ACs formed via both single 
and binary star channels. 
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1. Introduction 

Anomalous Cepheids (ACs) are short-period variable stars (from 
a few hours to two days), located on the colour-magnitude dia- 
gram ~ 1 mag above the horizontal branch, at colours similar 
to those of RR Lyrae stars ( Bono et al. 1997) . Their properties 
can be explained by assuming that they are relatively massive 
(1 -2M0), core-He burning stars in a very low-metaUicity regime 
that ignite heh um under partial electron-degeneracy conditions^ 



(Z < 0.0008, 'Renzini et al. 1977; Castellani & degl'Innocenti 
[Bono et al.. 1997: .Caputo et al.. ,2004: .Fiorentino et al 



1995 



2006) 



ACs have been observed in many nearby Local Group dwarf 
galaxies that were surveyed for short time variability, indepen- 
dently of the morphological type. Typically, nearby dwarfs host 
a few of them. They have been observed in purely old sys- 
tems {age > 9 - lOGyr), such as i n the dwarf spheroidal 
(dSph) s atellites Sculptor (3 ACs, 'Smith & Stryked Il986t 
Kaluznv et al. 1995), Sextans (6, Mateo et al. 19 95i), LeoII (4 , 



^iegel&Maiewski 2000), Ursa Minor (7. iNemec et all 1 19881). 
Draco llO. iHarris et al..l998 ':'Kinem uchi et al.l2 008'). an d also in 



isolat ed dSphs, like C etus (8 , Bernar d et al.ll200 9. Mone lli et all 
I2OI2I) . and Tucana (6, [Bernard et al.ll2009h . However, they have 



* Based on observations released by the Optical Gravitational 
Lensing Experiment (OGLE-III). 



also been found both in dSph galaxies with a large intermediate- 
age population (herei nafter, 1 - 6 Gyr), such as Fornax (17, 
lBersier& Wood* 2002), Cai'ina (15, lOalFOraet all 12003 1). and 
Leol (15, Hodge & Wright 19 78i) an d in gas-rich dwarfs such 
as Phoenix (12. iGallart et al.ll2004bh . The cases of Leo A and 
NGC 6822 are not as clear, f or the coexistence of ACs and short- 
period classical Cepheids dHoessel et al.l 119941: iDolphin et al.l 
I2OO2I : iBaldacci et ^1200 5). Furthermore, some have been iden- 
tified in four satelUtes of M31 ( Pritzl etal...2002 . 2Q04, 20051) 
for which the star formation history and the the population con- 
tent is still uncertain. Recently, three ACs have been propose d 
as candidates in the ultrafaint dwarf CVenI dKuehn et al.ll2008l) . 
However, nearby galaxies, especially the spheroidal satellites 
of the Milky Way, typically occupy large areas of the sky, so 
that probably the actual census is far from complete. On the 
other hand, ACs are very rare in globular clusters: so far, only 
one candidate has been confirmed in the metal-poor ([Fe/H] ~ 
-2 dex) NGC 5466 (Zinn & Dahn 1976), and a few others have 
been suggested (Corwin et aL, 1999: .Arellano Ferro et al]|2008l : 
lKuehnetal.ll20lilK 

The largest sample of ACs in nearby systems (83) has been 
collected for the LMC in the framework of the OGLE-III project 
(Optical Gravitational Lensing Experiment, Udalski et al. 200^- 
In a series of papers devoted to the analysis of thousands of 
variable stars in the LMC (classical Cepheids, ISoszvriski et al.l 



1 



G. Fiorentino & M. Monelli: AGs in LMG 



2008b|: ACs an d population II Ceph eids, ISoszvnski et alJ 
2008c; RR Lvrae. Sosz vnski et al.ll2009at long-period variables, 
Soszvnski et al. 2009b), the OGLE-III experiment released the 
V, I photometry, pulsation properties, and well-sampled light 
curves for variables with periods from a few hours to hun- 
dreds of days. In this paper, we take advantage of this rich data 
base to investigate the nature of ACs within the global picture 
depicted by the star formation history and chemical evolution 
available nowadays (e . g. [Ca rrera et al. 2008; Harris & Zaritskv 
l2009HSaha et al.l2010t ICioni et al.i201 luCarrera et al.i201 L) for 



the LMC. 

In fact, even though ACs identify a population of stars 
in a well-defined mass range, the origin of their progen- 
itors is still being debated. If they are the result of the 
evolution of single stars, then they trace a relatively young 
event in t he star formation, which occurred from ~ 1 to 
6Gyr ago dD emarque & Hirshfeld"1975'; 'Norris & Zinn) [19751 
ICastellani & degrinnocenti 1995; Caputo et al. 1999) . On the 
other hand, it has been proposed that they are the evolution of 
blue straggler stars (BSS), formed via mass-exchange in binary 
systems (Renzini et al. 1977; Sills et al. 2009) that survived in 
low-density environments, so ACs are tracers of the old, metal- 
poor population. The discovery of 83 ACs in a relatively metal- 
rich (Z - 0.008) environment such as that of the LMC makes 
this picture even more puzzling. 

In this work, we present the approach base d on pulsation 
mode l s built ad hoc for thi s class of variable stars ( Marconi et al.l 
120041; ICaputoet al] 12004 iFiorentino etal] 120061) to constrain 
the pulsation mass of individual stars. This, in turn, can be 
related to the star formation history (SFH) of the host sys- 
tem, as already succes sfully done to all the ACs discovered so 
far in nearby galaxies dCaputo et all 12004 iMonelli etani2012L 
IFiorentino e"tal]|2012l) . 

This paper is organized as follows. In Sect. |2] we give a 
full description of the sample selected from the OGLE-III sur- 
vey of the LMC and of the Fourier analysis needed to measure 
the pulsation amplitude in V-band. In Sect. [3] we compare the 
OGL E-III data with the th eoretical boundaries of the instability 
strip dMarconi et al. I I2OOI . In Sect.S we describe the theoreti- 
cal method used to simultaneously constrain the pulsation modes 
and the masses. Section |5]presents the radial distribution of dif- 
ferent samples of variable stars in the LMC. The discussion and 
the conclusions close the paper. 



2. The OGLE-III sample 

We selected the 83 AC s identified in th e short period Cepheids 
catalogue presented in ISoszvnski et alj (|2008c). The classifica- 
tion of the variable type is mostly based on the period-luminosity 
relation (PL), as detailed in Soszvnski et al. (2008a c). In partic- 
ular, the WesenheiQ index is a good diagnostic, because it in- 
cludes the colour information (and therefore the temperature) 
of the star without being affected by the reddening evaluation. 
Using this plane, it is straightforward for separating Cepheids 
(classical and anomalous) from population II Cepheids (here- 
inafter P2C, which include BL Her, W Vir, and RV Tau), be- 
cause th e latter have, at a given period, significantly fainter mag- 
nitudes (IFiorentino et al.l2006l and references therein). From the 
Wesenheit plane, it clearly emerges that and classical Cepheids 
occupy different regions. In particular, ACs are located at an 
intermediate position between classical Cepheids and P2C. 
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Fig. 1. Wesenheit (V, V - P) plane. Dots and triangles represent 
fundamental and first-overtone anomalous Cepheids as defined 
by OGLE-III (see Soszvnski et al. 2008c). The star represents 
one P2C candidate from OGLE-III that we have included in our 
sample. The solid line r epresents the theoret i cal pr ediction for 
BL Her type stars from iDi Criscienzo et a D (I2007h . while the 
dashed line is the observed Wesenheit relation for fundamental- 
mode classical Cepheids (.Soszvnski et al.. 2008 d) . 



However, the unambiguous distinction between anomalous and 
classical Cepheids is hampered by the fact that the first-overtone 
(FO) AC fall on the extension of the period-luminosity rela- 
tion to shorter p eriods of the fundamental ( F) classic al Cepheids 
(period < 2.5 d ICaputo et ai.ii2004; Soszvhski et aljr2008c ). For 
these stars, the shape of the Ught curve is the most efficient diag- 
nostic to separate the two groups. This picture is even more com- 
plicated by certain degeneracy existing with the RR Lyrae stars 
as well. Therefore, while ACs with periods longer than 1 d can be 
safely identified, in the short-period regime (< 1 d), brighter RR 
Lyrae stars and ACs can overlap on the period-magnitude plane. 
Moreover, the intrinsic similarities between the light curves of 
these two groups make the distinction uncertain. 

The final classification from OGLE-III includes 83 ACs. On 
the basis of their location in the Wesenheit plane, 62 were clas- 
sified as F and 21 as FO pulsators. We discuss the importance of 
the mode classification more fully in Sect. |4l To this sample, we 
added one star classified by OGLE-III as P2C (namely V-166), 
because it shares the pulsation properties of ACs, as we discuss 
in the next section. Figure [T] shows the (V, V - I) Wesenheit 
relation for the 84 stars. For comparison, w e also show the the- 
oretical relation for BL Her stars given bv iDi Criscienzo et al.l 
( 2007) . which are the subsample of the P2C class with a short 
period (P < 4 d). The Wesenheit relati on for fundamen t al-mod e 
classical Cepheids has been taken from lSoszvnski et alJ(l2008al) . 

In Table [1] we summarize some of the main properties pro- 
vided by the OGLE-III catalogue: the first eight columns give, 
for each variable star, the identification number, the position 
(right ascension and declination), the mean V and / magnitudes, 
the period, the amplitude in the / band (A/), and the mode clas- 
sification. The V-band amplitude. Ay, was not given because the 
number of phase points in the V-band (27 to 138) is typically a 
factor ~ 10 smaller than in / (319 to 1299). Nevertheless, the 
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minimum number of data points in the V band is 27 epochs, en- 
suring well-sampled light curves. For this reason, we performed 
the Fourier analysis on V-band photometry to give an indepen- 
dent estimation of the pulsation parameter and provide the am- 
plitude A i/. Taking advantage of the accurate periods given from 
OGLE-III catalogue, we estimat ed Ay using the GRAT IS code 
(developed by P. Montegriffo, see lFiorentino et al.l2010i and ref- 
erence therein). The results are reported in columns 9 and 10 
of Table [1] The comparison with the OGLE-III mean magni- 
tude reveals that this difference A < V > is always less than 
about 0.05 mag. The total error budget on our evaluation gives 
o-(Av) S 0.1 mag, thus perfectly suitable for the goals of the 
present analysis. The star V-82 was excluded because no V-band 
photometry is available. 

3. Comparison with theoretical predictions 

To properly compare theory and observations, we need to ac- 
count for the distance modulus and the reddening of the host 
galaxy. The distance modulus to the LMC has been exten- 
sively discussed in the literature, since this galaxy is the cor- 
nerstone of the distance scale. We adopt here fiQ = 18.50 mag, 
as given by a recent estima te using Red Clu mp stars observed 
in optical and near-IR (Fio rentino et alj|201 Ih . This agrees with 
many of the measurements d erived by independent methods, in- 
cluding classical Cepheids ('Bono el al. 2002), RR Lyrae stars 
(Catelan & Cortes 2008), planetary nebulae luminosity function 
tReid & Parker 2010), the red clum p, and the tip of the red gi- 
ant branch (Romaniell o et al.ll2000h . The AC sample distributed 
across the whole area surveyed by OGLE-III, we cannot adopt a 
unique reddening value for all the variables. To correct the mag- 
nitude of each star for the appropriate extinction value, we u sed 
the reddening map recently derived bv lHaschke et al.l ( 1201 lb us- 
ing red clump and RR Lyrae stars from the same OGLE-III re- 
lease. In particular, we assumed the reddening value, E{V - I), 
for the closest position available for each individual AC. This 
is reported in column 11 of Table [1] The mean reddening is 
E{V - I) ~ 0.08 mag, and only one star has a value higher than 
0.2 mag. 

In Fig. |2]we show the comparison between the theoretical 
prediction for the bou ndaries of the pulsation instability strip 
dFiorentino et alj|2006h and the sample of ACs. The top and bot- 
tom panels show, as a function of the logarithm of the period, 
the absolute mean colour (< My > - < M/ >) and the ab- 
solute mean < My > magnitude, respectively. The agreement 
between theory and observations is good in both cases and for 
all the objects (the symbols are the same as in Fig.lTJ. The star, 
classified by OGLE-III as population II Cepheids, seems to fully 
agree with AC observational properties. We note that there are 
few outliers in the upper panel. This is possibly due to the uncer- 
tainty in the reddening determination, with directly affects the 
true colour of stars. Nevertheless, this comparison suggests that 
the correction for reddening applied is satisfactory for most of 
the objects. 

4. Constraining masses and pulsation modes 

AC stars follow well-defined relations, which are the period- 
mass-amplitude relation (PMA) and the mass dependent period- 
luminosity-colour relation (hereinafter PMC). Here, it is worth 
mentioning that this theoretical scenario was used for classical 
Cepheids to constrain pulsation masses with very high preci- 
sion (1%) and r eturned values in good agreement with dynami- 
cal masses (e.g. lPietrzvnski et al.ll2010t rPietrzvnski et al.il201 ll 
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Fig. 2. Comparison between theoretical models and OGLE- 
III observations. The distance modulu s assumed for LMC is 
= 18.50 mag dFiorentino et al.ll201 ih . and the individual red- 
den ing values are as sumed using the reddening map given by 
Haschke et al.l (1201 lb (see text for details). The coding of the 
colours is the same as in previous figures. Top - (< My > - < 
Ml >)() colour vs period (log P) diagram are compared with the- 
oretical blue (solid line) and red (dashed line) edges of the insta- 
bility strip (Marconi et al.ll2004l) : Bottom - Absolute magnitude 
(My) VS period (log P) plane compared with the theoretical first 
overtone blue (solid line) and the fundamental red (dashed line) 
edges of the instability strip. 
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Fig. 3. Mass values for anomalous Cepheids in the LMC as de- 
rived from PMA and PMC relationships. Solid grey line shows 
the relation of equality. The coding of the colours is the same 
as in previous figures. The ACs identified as FO in OGLE but 
F in this work are highlighted with red open circles, while the 
ones classified as F in OGLE and FO in this work have been 
emphasized with open squares. 
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Table 1. The anomalous Cepheids from the OGLE-III sample, including 1 star (last row) classified by OGLE-III as P2C (see text 
for details.) 



OGLE-in 



This work 



ID 



R.A. 

J2000 



Decl. 
J2000 



V 

mag 



I 

mag 



A, 
mag 



modeoGLE 



<V> 
mag 



Ay 

mag 



EiV-I) 
mag 



Mf 
Mo 



Mro 
Mo 



modem 



001 

002 
003 
004 
005 
006 
007 
008 
009 
010 

on 

012 
013 
014 
015 
016 
017 
018 
019 
020 
021 
022 
023 
024 
025 
026 
027 
028 
029 
030 
031 
032 
033 
034 
035 
036 
037 
038 
039 
040 
041 
042 
043 
044 
045 
046 
047 
048 
049 
050 
051 
052 
053 
054 
055 
056 
057 
058 
059 
060 
061 
062 
063 
064 
065 
066 
067 
068 
069 
070 
071 
072 
073 
074 
075 
076 
077 
078 
079 
080 
081 
082 
083 
166 



04:37:18.96 
04:41:56.22 
04:43:16.03 
04:47:21.63 
04:51:11.62 
04:52:34.51 
04:57:31.48 
04:58:24.59 
04:58:51.29 
04:59:00.10 
04:59:38.09 
04:59:49.43 
04:59:56.62 
05:00:08.26 
05:00:59.32 
05:01:36.69 
05:02:03.13 
05:02:36.36 
05:03:17.71 
05:05:54.54 
05:06:37.49 
05:06:45.08 
05:07:52.42 
05:08:44.05 
05:09:55.67 
05:10:42.62 
05:10:45.77 
05:13:10.02 
05:13:33.58 
05:15:07.61 
05:15:32.98 
05:15:56.14 
05:16:08.01 
05:17:11.15 
05:18:22.19 
05:18:58.84 
05:19:16.66 
05:19:59.97 
05:20:44.46 
05:21:13.22 
05:21:14.35 
05:23:34.60 
05:24:35.40 
05:25:54.11 
05:26:24.68 
05:26:27.17 
05:27:05.27 
05:27:12.12 
05:28:03.57 
05:28:57.71 
05:30:14.24 
05:31:01.53 
05:31:06.20 
05:31:06.72 
05:31:41.11 
05:31:49.45 
05:31:49.88 
05:33:39.38 
05:34:14.81 
05:36:08.89 
05:36:24.32 
05:36:58.43 
05:37:54.39 
05:39:26.03 
05:40:03.04 
05:45:53.29 
05:48:22.07 
05:51:09.02 
05:53:31.91 
05:53:43.23 
05:54:43.24 
05:59:08.23 
06:00:04.99 
06:02:45.18 
06:03:33.33 
06:04:33.02 
06:04:35.73 
06:06:58.20 
06:07:02.01 
06:09:35.12 
06:09:38.40 
06:16:58.31 
06:18:52.45 
05:16:58.99 



-69:49:11.5 
-66:51:59.1 
-69:07:57.6 
-70:22:12.0 
-69:00:33.0 
-69:43:27.5 
-70:15:53.6 
-71:05:11.5 
-67:44:23.9 
-68:14:01.1 
-70:37:45.5 
-69:33:12.4 
-70:42:24.9 
-67:54:04.3 
-69:33:24.9 
-67:51:33.1 
-68:09:30.4 
-69:40:12.2 
-68:49:33.5 
-70:43:30.3 
-68:23:40.3 
-70:52:44.7 
-68:50:28.8 
-68:46:01.2 
-70:03:08.2 
-68:48:19.6 
-67:10:05.6 
-68:46:11.9 
-67:33:42.2 
-68:04:21.4 
-68:10:29.4 
-69:01:29.1 
-70:39:16.5 
-69:58:33.1 
-69:03:38.4 
-69:26:47.8 
-70:11:58.4 
-71:45:23.7 
-69:47:46.5 
-70:34:20.7 
-70:29:39.5 
-69:10:58.2 
-68:48:22.6 
-69:26:52.9 
-68:57:55.0 
-69:58:57.0 
-71:23:33.4 
-69:37:19.6 
-69:39:15.2 
-70:07:15.5 
-68:42:31.1 
-70:42:22.2 
-68:43:45.3 
-68:22:29.8 
-68:44:37.7 
-70:33:22.6 
-70:46:30.0 
-70:21:28.6 
-68:23:43.0 
-70:37:01.3 
-71:30:39.4 
-70:46:08.1 
-69:19:28.7 
-71:24:47.6 
-70:04:47.8 
-72:01:19.6 
-70:45:49.3 
-70:45:42.6 
-71:32:44.8 
-69:17:06.0 
-70:10:16.1 
-68:24:01.2 
-70:44:59.8 
-70:45:31.2 
-68:27:41.4 
-70:57:55.8 
-71:40:35.8 
-72:52:08.7 
-69:31:55.2 
-70:10:42.1 
-69:34:04.1 
-70:52:18.9 
-70:51:55.4 
-69:51:19.3 



18.676 
18.268 
18.794 
17.719 
18.526 
17.671 
18.247 
17.911 
17.913 
18.723 
18.254 
18.663 
18.488 
17.241 
17.587 
17.448 
18.190 
18.213 
17.891 
18.359 
17.827 
18.378 
17.763 
17.986 
18.413 
17.483 
17.558 
18.322 
18.745 
17.685 
17.727 
17.780 
17.308 
18.474 
18.658 
17.787 
17.743 
17.908 
18.245 
18.037 
18.200 
18.715 
18.593 
17.609 
18.976 
17.851 
17.482 
17.324 
18.451 
17.049 
18.845 
17.577 
17.299 
18.798 
17.603 
17.877 
17.455 
18.055 
20.004 
17.527 
18.393 
17.659 
18.756 
17.619 
17.508 
18.233 
18.451 
18.809 
17.460 
18.370 
17.755 
17.535 
17.396 
17.551 
18.490 
17.671 
18.099 
17.417 
17.634 
17.997 
18.504 

18.236 
17.696 



18.018 
17.621 
18.313 
17.018 
17.858 
17.057 
17.688 
17.300 
17.344 
18.068 
17.671 
17.969 
17.963 
16.639 
16.962 
16.926 
17.585 
17.562 
17.414 
17.998 
17.188 
17.796 
17.194 
17.679 
17.904 
16.816 
16.956 
17.574 
18.120 
17.258 
17.215 
17.167 
16.616 
17.874 
18.143 
17.160 
17.132 
17.294 
17.660 
17.433 
17.625 
17.897 
17.913 
17.052 
18.325 
17.264 
16.881 
16.718 
18.011 
16.609 
18.200 
17.008 
16.738 
17.901 
17.011 
17.284 
16.813 
17.710 
18.724 
16.963 
17.745 
17.162 
18.013 
17.011 
17.041 
17.585 
17.786 
18.265 
16.823 
17.735 
17.330 
17.038 
16.781 
16.889 
17.876 
17.032 
17.459 
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Fig. 4. Top- Same as Fig. [T] but we highlight the stars for which 
we obtain a different mode classification as in Fig. [3] Bottom - 
Amplitude {Ay) vs period (logP) diagram. The coding of the 
colours is the same used in Fig. [T] The amplitudes have been 
derived in this paper 

ICassisi & SalarislbOllh . In particular, it is possible to estimate 
individual masses once distance and reddening are kn own by 
applying the following equations, adapted from Marconi et al.l 
(l2004 : 

Mf^pmaIM^ = 10(0.013-0.53 <M.>-1.3 logF-0.244 A,) 
Mf,PMcIM^ = 10(-° '''-°-5^ <M,>-I.55 logF+I.44«Mv>-<M,») 
Mfo.PMc/M^ = 10(-1.29-0.58 <M„>-1.79 log/'+1.54«M„>-<M,») 

where My, Mj are the absolute intensity-weighted mean V, 
I magnitudes, P is the period, and Ay is the amplitude in 
the V-band. In this work, we assumed reddening law from 
ICardelli et alj (Il989h . The formal uncertainty on the derived 
masses is 5% usin g the PMC and 15% using the PMA 
(iMarconi et al ] 12001 . It is worth mentioning that the compar- 
ison of the masses derived with a different relation is indepen- 
dent of the distance modulus adopted, but it does depend on the 
reddening correction applied. Assuming the OGLE-III classifi- 
cation, we estimated the mass of each star, which is reported in 
Table[T] We provide the mean mass from Eqs. 1 and 2 in the case 
of F, and the mass from Eq. 3 for FO. For three objects, we were 
not able to estimate a reliable mass, e.g. the mass estimated for 
the FU V-24 and V-59 with Eqs. 1 and 2 shows high discrepan- 
cies, and the mass estimated for FO V-28 is too high (~ 2.3M0) 
for an AC. 

A method is discussed in iMarconi et a n (12004 which uses 
the previous relations to simultaneously constrain the pulsation 
mode and the mass of the star This method is based on the fact 
that the PMA relation is only valid for the F pulsator, whereas 
the PMC exists for both pulsation modes. Therefore, whenever 
the PMA and the PMC for F pulsators give consistent mass, 
this is a robust indication that the star is actually pulsating in 
the fundamental mode. The correct mode classification is im- 
portant because it affects the mass determination. For the same 
luminosity and colour, we expect F pulsator to be more massive 



than FO pulsators. However, the usual diagnostics used for other 
classes of variable stars, such as the Bailey diagram (Ay vs log P 
plane) or the light curve morphology, do not allow the pulsation 
mode to be uniquel y constrained. For th is reason, we applied 
the method given in IMarconi et al.l (l200l to check whether the 
mode classification based on the Wesenheit index (provided by 
the OGLE team) is consistent with theoretical pulsation predic- 
tions. 

Figure |3] compares the mass values obtained using Eqs. 1 
and 2. The grey line shows the relation of equality. Assuming 
the formal error for the theoretical relations, we classified as F 
pulsators all the stars for which the mass difference is within 
the sum of the two 1 - cr errors. For a better identification, we 
have also classified F pulsators all those stars that did not pass 
the previous criterion, but show high PMC masses ( > 2.5Mq), 
which are unrealistic in the adopted theoretical scenario. Out of 
83 stars, we then classify 66 F and 17 FO, as reported in the last 
column of Table[T] Three stars (V-24, V-28, and V-59) present a 
large mass difference and also an atypical amplitude ratio, and 
at least in one case, a very red colour This suggests that they 
might be blended sources, so they are reported as unclassified 
in Table [T](n.c.). Moreover, V-166 previously identified as P2C 
turns out to be an AC pulsating in F mode. 

The comparison with the OGLE-III classification reveals 
some discrepancy. Out of 19 FO for which we could estimate 
the pulsation mode, we only recover 12. The remaining seven 
stars were classified as fundamental mode pulsators with our 
method, and are labeled as "F**" in Table [T] and highlighted in 
Fig. |3] For these stars, we also report, in parenthesis, the mass 
value corresponding to the F mode given by our classification. 
For these stars the FO mass value has a reliable value, systemat- 
ically > IMq. Then, masses of these "F**" stars are compatible 
with both F pulsators with relatively higher mass, and FO pul- 
sator of lower mass. 

On the other hand, there are three stars, flagged with "FO*" 
in Table [T] and classified as FO, differently from the OGLE-III 
work. For these stars, the masses derived with the PMC re- 
sult to be always higher than the ones derived with the PMA 
less than 1 ~ Miti)^ in disagreement with the general behaviour 
Caputo et al.ll2004]) . Their position in Fig.|3]shows that they are 
close to the bulk of F-mode stars. 

To investigate the reason for these differences, we recalcu- 
lated the masses of the 12 outliers varying the values of the red- 
dening and the amplitude Ay (which could be affected by blend- 
ing). We found that, while modifying the amplitude has little or 
no effect, modifying the reddening value has a stronger effect. As 
an example, increasing the reddening of 0. 1 changes the classifi- 
cation of 2 FO (FO* in table) out of 3 to F and 5 F (F** in table) 
out of 7 to FO. This is an indication that the uncertainty on the 
reddening is the dominant source of error in the present analysis. 

In Fig. |4] we show the Wesenheit relation, similar to Fig. [T] 
(top) and the P - Ay (bottom), and we highlight in both panels 
the stars for which we find a different mode classification, as in 
Fig.|3] We note that the FO our method recovers as F all lie on the 
lower envelope of the FO relation in the Wesenheit plane, thus 
supporting the classification's consistency with both FO and F 
modes. Also, the Bailey diagram clearly shows that, as opposed 
to the case of RR Lyrae stars, this plane is not an efficient diag- 
nostic to separate F and FO pulsators. In fact, while the shortest 
period ACs are clearly FO, while those with large amplitudes 
tend to be F, there is a clear degeneracy in the low-amplitude 
regime, for periods with -0.2 < logP < 0.1 d. 
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Fig. 5. Spatial distribution of anomalous (red), population II (blue), and classical Cepheids (black), RR Lyrae (grey) in the total area 
of the Large Magellanic Cloud (~ 39.7 square degrees) covered by OGLE III survey. 
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Fig. 6. Cumulative distributions for anomalous (red), population 
II (blue), and classical Cepheids (black), RR Lyrae (grey). 

5. Spatial distributions of variable stars 

To better understand the origin of ACs we take advantage of the 
extensive and complete sample of variable stars discovered by 
OGLE-III survey of the LMC, namely 22,651 RR Lyrae {ab and 
c type), 203 P2C, and 3,087 young classical Cepheids (F and 
FO only). The Magellanic clouds are the only galaxies where 
such an extensive variable star study has been performed, cov- 
ering variables with a wide range of periods. The total area cov- 
ered of the LMC by the OGLE-III project is 39.7 square degrees 
and samples the central regions of the galaxy. In particular, it 
covers the bar and the inner disc out to a maximum distance 



of ~ 5 kpc. However, the distribution of OGLE-III fields fol- 
lows the bar direction from north-west to south-east. Also we 
note that the size of the LMC is significantly g reater than the 
area covered by OGLE-III (e.g. ' Saha et aD20ltf ). Nevertheless, 
the present data set allows studying the spatial distribution of 
the different samples of variable stars over the central region 
of the galaxy, thus tracing the parent population. As a matter 
of fact, the RR Lyrae stars and the P2Cs are horizontal branch 
stars tracing the old and predominantly metal-poor population 
([Fe/H] < - 1, dex), while the classical Cepheids trace the 
young and mainly metal-rich ([Fe/H] > - 0.5 dex) stars. As 
shown in Fig. |5] the spatial distribution of young stars is con- 
centrated at the centre of the galaxy, nicely tracing the bar re- 
gion where the star formation occurred at very recent epochs 
( < 0.5 Gyr ago). Also, there are two clear overdensities at the 
end of the bar at ~ 85 and 74 degrees. On the other hand, the RR 
Lyrae stars seem to follow a much broader uniform distribution 
as is usually observed in other galaxies (Monelli et al. 2010a,^ 
l2012Hde Boer et al.ll201 ll) . Interestingly, the ACs do not seem to 
obviously follow either of the two previous distributions. In par- 
ticular, while the majority seems concentrated in the innermost 
regions, where the bulk of the RR Lyrae stars is located, they are 
not usually placed in the bar Moreover, there is a clear overabun- 
dance in the east direction (RA > 87 degrees), where the density 
of both RR Lyrae and Cepheids has significantly decreased. A 
similar excess does not seem to be present on the other side of 
the galaxy. 

To better compare the spatial distribution of the different 
class of variable stars, we built the cumulative distributions as a 
function of the gala ctocentric radius using the LMC model from 
Ivan der Marell (1200 L) to derive the deprojected distance from the 
centre. They are shown in Fig.|6]for the RR Lyrae stars, classical, 
P2Cs, and ACs. We have also computed, using a Kolmogorov- 
Smirnov test, the probabilities that the various sample are drawn 
from the same parent population as the RR Lyrae stars. We found 
that this is 0% (5 x 10"''%) for classical Cepheids, 0.2% for ACs, 
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Fig. 7. AC frequency (log 5) per 10^ Lq as a function of the ab- 
solute V magnitude of the host galaxy (My). Open circles show 
purely old dSph galaxies, while starred symbols indicate galax- 
ies with large intermediate-age populations. The crosses mark 
the four M31 satellites. The LMC and SMC are marked by the 
open triangle (for an assumed luminosity of 2 x lO^L©) and 
square, respectively. The solid line show the fit to the relation for 
the old systems only, while the dashed line is the fit to the four 
crosses. We note that the present value for the LMC is a lower 
limit, since ACs have so far been detected only in the central 
region of the LMC. However, the LMC is known to be signifi- 
cantly larger, possibly occupying an area ~ 9 times larger (see 
e.g. ,Irwia.l99k.Maiewski et al...2009;.Saha et al...201Q) . 

and 87% for P2C. Therefore, the P2Cs seem to follow the same 
distribution of the RR Lyrae stars, in agreement with the fact that 
they both originate in the old population. We have also com- 
puted the probability that the ACs come from young classical 
Cepheids, which is 8%. To avoid possible bias from the area 
covered by the OGLE-III survey, we repeated these estimations 
using only subgroups of stars located at distances within 2, 3, 
and 4 kpc from the LMC centre. We find that the probability that 
the distributions of the RR Lyrae and P2C coincide is always 
above 75%, while that of RR Lyrae and classical Cepheids is 
always below 1%. In the case of ACs and RR Lyrae stars, the 
probability widely fluctuates but, with some hint of decreasing 
probability for increasing the limit radius, ranging from 32, to 
10, to 0.5% in the three cases. 

These tests do not show any clear similarity between the dis- 
tributions of ACs with either that of classical Cepheids or RR 
Lyrae stars. This suggests that they are not uniquely correlated 
either with the old or with the youngest LMC population. 

6. Discussion of the origin of ACs 

The theoretical scenario that explains the structure, evolution, 
and pulsation properties is settled for different classes of pul- 
sating variable stars. This has strong implications in the study 
of the star formation history of galaxies, because it is possible 
to characterize the properties of the parent population (e.g., age 
and metallicity) by investigating the bright variables associated. 
However, the presence of ACs does not uniquely identify the 



nature of their progenitors. The two most favoured formation 
channels are the evolution of single, metal-poor star with mass 
< 2Mq, or the evolution of coalescent binary systems of metal- 
poor stars. Independent of the formation mechanism, a crucial 
characteristic of the progenitor for producing an AC appears to 
be its low metallicity. In fact, central-helium burning structures 
in the AC mass range, but more metal-rich than Z ~ 0.0008 
dpallart et al. 2004a; Fiorentino et al. 2006), are too red to en- 
ter the Cepheid instability strip and populate the red clump. To 
have pulsators at higher metallicity, the mass of the progenitor 
must be higher. In this case, the star will cross the strip at higher 
luminosity producing short period classical Cepheids, which ig- 
nite the helium quiescently. In particular, the minimum mass of 
central helium-burning pulsators crossing the instability strip in- 
creases fr om L9Mo (for Z = 0.0004) to 3.6Mq (for Z = 0.008, 
ICaputo e t al. 2004). 

The binary origin is a plausible explanation of ACs having 
been ubiquitously discovered in low-density, very metal-poor 
environments (Z < 0.0008), such as nearby purely old dwarf 
galaxies. In fact, ACs have been observed in dwarf galaxies that 
do not show any intermediate-age population, such as Draco, 
Ursa Minor, Sculptor, Sextans, Tucana, LeoII, and Cetus. This 
implies that, at least in these galaxies, ACs must be the progeny 
of metal-poor, binary systems. On the other hand, very few sys- 
tems have been discovered to host both classical and anoma- 
lous Cepheids. Gallart et al. (2004b) first pointed out this oc- 
currence in the case of the transition dIrr/dSph galaxy Phoenix, 
supporting the evidence that the same applies to Sextans A and 
Leo A. Comparing the spatial d i stribu tion of different samples 
of stars in Phoenix. iGallart et al.l te004b) find that, while classi- 
cal Cepheids are strongly segregated in the central regions, fol- 
lowing the spatial distribution of the youngest population, the 
ACs are distributed over a broader area, but are still more con- 
centrated than the RR Lyrae stars, which can be safely taken 
as representative of the old population. This suggests that the 
ACs in Phoenix are mostly intermediate-age stars, though the 
possibility that a few are BSS descendents cannot be excluded. 
This implies that, despite its prolonged star formation until 
epochs that are recent enough to produce classical Cepheids, 
Pho enix experienced lit tle chemical enrichment, in agreement 
with lHidalgo et al.l (12009:) . and the 1-6 Gyr old population was 
still metal-poor enough to produce ACs. 

6.1. The case of LMC 

The case of the LMC presented here is probably more complex, 
since the LMC is two orders of magnitude more massive than the 
typical satellite dSph, and has a comp lex structure (i.e. disc an d 
bar). Nevertheless, until the work by ISoszvhski et al.l (l2008cl) . 
no ACs had been unambiguously detected. The present sam- 
ple of 84 objects is the largest so far for any external galaxy. 
The coexistence with large samples of different variable stars al- 
lowed us to compare the radial distribution of these objects using 
a Kolmogorov-Smirnov test. We found that the present sample 
of ACs cannot be associated to the young classical Cepheids. 
It must be stressed here that given the large sample of AC and 
classical Cepheids in the LMC, their different nature clearly ap- 
pears for the first time. Despite ACs occupy the sa me instabil- 
ity stri p of the classical Cepheids as suggested by iCaputo et al.l 
(|2004|) and populate the low-mass and low-metallicity exten- 
sion, it is clear that the two groups follow different Wesenheit 
relations (Soszvhski et al. 2008c). We also found that the radial 
distribution of ACs do not correlate with the oldest population, 
as represented by RR Lyrae and P2C, suggesting that, as for of 
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case of Phoenix (iGallart et al.ll2004bl) . a large fraction of single 
intermediate-age stars is present. 

To investigate this hypothesis, we interpreted the properties 
of ACs in the full picture of the LMC star formation history. 
Given the mass derived in this work (see Sect. |4]i, which is 
1.2 + Q.2Mq for the bulk of ACs, we can estimate the corre- 
spondin g age. Using scaled sol ar evolutionary tracks from the 
BaSTfll Pietrinferni et al.l20 04') database, we find ages from 1.8 
to 5.8 Gyr, assuming Z=0.0006. However, the age-metallicity 
relation (AMR) ^ based on spect r oscopic analysis, as prese nted 
in recent years dCole et al.ll2005l : ICarrera et al. l l2008ll20Tll) . re- 
veals that the mean metalUcity in this age range is an order of 
magnitude greater than the expected maximum value for ACs, 
being close to Z-0 .008. Interestingly , the s patial analysis by 
ICarrera et all (|2008|) and ICarrera et all (1201 ll) reveals that, in a 
number of disc fields at diff'erent distances from the LMC cen- 
tre, but also in the bar field from Cole et al. (2005), the metal- 
Ucity distribution always presents a well-populated tail toward 
low metallicity (Z ~ 0.0004), but the most metal-poor stars 
are systematically older than 10 Gyr. This occurrence seems at 
odds with our hypothesis that the bulk of ACs are metal-poor 
intermediate-age stars. However, one could ask whether, due to 
the limited number of spectroscopic targets so far analysed, we 
are facing an observational bias, and that very few metal-poor 
stars exist in the age range from ~ 2 - 6 Gyr. 

We can attempt to estimate the percentage of the mass in- 
volved in the star formation episode that generated the ACs, fol- 
lowing the simple prescriptions given in Renzini & Fusi Pecc^ 
(1 19881) . which hold for each post-main sequence phase. The 
number of ACs (Nac - 84) observed is related to the total lu- 
minosity (Lz<o.ooo8) of the population they belong to via the fol- 
lowing relation: Nac - B(t)Lz<uMm'''- Where B(t) is the spe- 
cific evolutionary flu x (~ 0.15 x lO"'" stars per Lojr"', see 
iRenzini & Fusi Peccii ri988). t is the time spent in the instability 
strip f or a mass of about 1.2Mq (~ 50 Myr, see iFiorentino et aU 
I2006h . andLz<o.ooo8 results to be IxIO^Lq- This value is very low 
related with the total luminosity of the LMC (Ltot - 2x IO^Lq), 
even if we account for the limited area covered by the OGLE- 
III survey. Therefore, the corresponding expected total number 
of stars from this population is very low, a few thousand over 
the whole field-of-view, implying that few stars brighter than the 
Horizonal Branch are expected per square degree. This strongly 
suggests that such a small event cannot be properly sampled by 
the SFH analysis or spectroscopic investigation of incomplete 
samples. This in turn highlights the key role of ACs in this con- 
text. In fact, the discovery that a small amount of gas was still 
poorly enriched 6 Gyr ago gives additional constraints to the star 
formation history that have to be accounted for in the formation 
and evolution modeling of the LMC. 

Finally, we want to briefly comment the possibility that 
ACs are young, but metal-rich stars looks unfeasible. Stars with 
metallicity close to Z = 0.008 are expected to be significantly 
more massive to enter the instability strip during the central 
helium-burning phase, with respect of more metal-poor stars. 
Even if we invoke unusually high mass loss to reconcile with the 
low-mass we estimated in this work, the mass of the helium core 
would not change significantly, leaving us with brighter stars, 
which would be at odds with the observed properties of ACs. 



^ http://albione.oa-teramo.inaf.it/ 



6.2. Comparison with otiier galaxies 

ACs appear to be intrinsically rare objects. This is valid in gen- 
eral for nearby dwarf galaxies and for the LMC in particular. The 
OGLE- III projects surveyed more than 35 million stars, and only 
84 ACs were detected. In this case, we can safely assume that the 
completeness is very high, both because of the large number of 
measurements and the photometric quality of the data. 

Nevertheless, Mateo et al.l (119951) first notice that a corre- 
lation exists between the frequency of ACs (log 5, per 10^ Lq) 
and the total luminosity of the host galaxy (My), valid for many 
nearby satellites. In particular, it was found that the frequency 
of ACs decreases for increasing luminosity of the host galaxy. 
Prit zl et al.l (|2002|) noticed that th e surveyed M3 1 sat ellites also 
follow this trend, as does Phoenix (iGallart et al.l2004"b ). Figure]?] 
shows the frequency of ACs as a function of the absolute vi- 
sual magni tude of the host galaxies. We used the data from 
Table 5 of iPritzl et al.l (2004), to which we added the values 
for Phoenix p allart et al. 2004 bl), Cefiis (Bernai'd etai]|2009l 
lMonellieta"l]|2012i) . Tucana (Bernard et al.i ,2009,) . Andromeda 

I, and Andromeda III (Pritzl etal. 2005). The symbols used 
in the figure diff'erentiate the objects according to their SFH 
(see the caption for details). In particular, we represent purely 
old systems (i.e. Ursa Minor, Draco, Sculptor, Leo I, Sextans, 
Tucana, and Cetus) and dwarf galaxies characterized by an im- 
portant intermediate-age population (Carina, Fornax, Phoenix, 
and LeoII). We also plot four M31 sateflites (And I, And II, 
And III, and And VI), for which the data available are not deep 
enough to set tight constraints on their SFH. In fact, the pho- 
tometry reaches roughly one mag below the horizontal branch, 
and it is impossible to clearly detect the presence of a strong 
intermedia te-age population (fo r example, compare with Carina 
in Fig. 1 of lHarbeck etani2001h . 

The two lines are the fits to the two group of galaxies. The 
plot suggests that purely old systems and intermediate-age pop- 
ulation systems follow different relations. In particular, when 
an intermediate-age population is present, the frequency of ACs 
tends to be higher than expected from a purely old population, 
and the discrepancy increases for increasing luminosity. The 
underlying physical mechanism of such a relation is unclear. 
However we only expect ACs resulting from the evolution of 
binary systems in purely old systems. Thus, the correlation with 
the luminosity, hence the mass, of the host galaxy may indicate 
that binary systems have higher chances of surviving in low- 
mass systems. It is worth recalling in this cont ext that a simi- 
lar coi Telation has been found for BSSs. In fact. iMomanv et al] 
(l2007l) have first shown that the frequency of BSS with respect 
of the number of horizontal branch stars decreases for increas- 
ing luminosity of the host galaxy. This might be an independent 
evidence that an evolutionary link exists between the binary sys- 
tems, BSSs, and ACs. If such a correlation could be explained 
this way for purely old systems, the fact that galaxies known 
to have intermediate-age population present an excess of ACs, 
compared to purely old galaxies, would be simply due to the co- 
existence of ACs formed via both formation channels. Among 
the four M3 1 satellites, two follow the fit of the old dSphs (And 

II, And III), while two present an apparent excess of ACs (And 
I, And VI). This is possibly because the small fraction of these 
galaxies covered by the data for these galaxies introduces some 
bias, but it is also possible that And I and And VI host some 
intermediate-age population, at present undetected but responsi- 
ble for this effect. 

What seems really surprising is that both Magellanic clouds 
appear to follow the same relation of these four galaxies, but in 
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a significantly higher luminosity regime. It is also interesting to 
note that the SMC hosts very few ACs candidates (3 F and 3 FO 
[Soszvnski et al. 2010). As for to the LMC, it presents an excess 
of ACs, thus supporting that it is possible that a small population 
of intermediate but metal-poor stars exists. 

The new OGLE-IV survey will possibly shed Ught on this 
problem. In fact, it is well known that a significant age gra dient is 
prese nt in the disc of the LMC (Gallart et al. 2008; Carre raet alJ 
1201 lb in the sense that the age of the youngest population gets 
older for increasing radius. This means that, starting from a cer- 
tain distance, we do not expect any population that is young 
enough to produce ACs, thus leaving only the BSS aftermath, 
if present. Moreover, the absence of young main sequence stars 
will permit direct detection of BSSs, thus allowing a direct cor- 
relation with ACs. The distance where this is expected to occur 
is far beyond the present coverage of the OGLE-III survey, but 
it will be definitely covered by OGLE-IV, at present in execu- 
tion. Stronger constraints on the formation mechanisms of the 
ACs in the LMC are therefore expected in the near future. Also, 
a systematic investigation of the M3 1 satellites to derive the de- 
tails of their SFH, which is at present perfectly feasible for the 
ACS and WFC3 cameras on-board the Hubble Space Telescope, 
would help solve this problem. 



7. Conclusions 

We have analysed the pulsation and spatial properties of 83 ACs 
discovered by the OGLE-III survey in the LMC. The main re- 
sults obtained are the following: 

• We classified one star (V-166), previously identified as 
P2Cs, as bona-fide ACs, thus increasing this total number to 84 
stars. 

• We performed the Fourier analysis of the V-band data, de- 
riving the V-band amplitude Ay and mean V values for the 83 
stars for which the V photometry is available, with only one 
excluded (V-82). W e then applied the method introduced by 
iMarconi et al.l (l2004l) to simultaneously derive individual masses 
and pulsation modes. We find good agreement with the OGLE- 
III classification for most of the stars (~ 85%), based on the 
Wesenheit relation and the morphology of light curves. The mass 
of the ACs ranges from 0.8 to 1.8 Mq, with a mean value of 
1.2 + 0.2Ma. 

• The large sample of ACs allowed, for the first time, to un- 
ambiguously show that anomalous and classical Cepheids are 
different. In fact, even though they populate the same instability 
strip, they follow distinct period-luminosity relations and differ- 
ent spatial distributions. 

• We compared the spatial distribution of the ACs with 
that of other samples of variables (RR Lyrae, P2C, classi- 
cal Cepheids) using a Kolmogorov-Smirnov statistical test. We 
found that RR Lyrae and P2C, both bona-fide representatives 
of the old population (> lOGyr) follow the same distribution, 
which is different from that of the young classical Cepheids. The 
distribution of ACs does not resemble any of the other samples. 
We conclude that a large fraction of the ACs descend from sin- 
gle, metal-poor, intermediate-age stars that are 1-6 Gyr old. 

• The comparison with the presently available SFH and 
AMR shows that such a population is negligible. We show 
that the total mass formed in that star formation event is about 
5 X lO"'' the total mass of the LMC, thus resulting in very few 
stars that are very unlikely to be observed in spectroscopic sur- 
veys or in statistical approaches to derive the SFH. This shows 
the importance of variable stars and rare objects like ACs for 



giving independent constraints on the modeling of galaxy evolu- 
tion. 

• We reanalysed the f requency-luminosity relation discov- 
ered bv lMateo et al.l (Il995h . showing that, if only purely old sys- 
tems are considered, the correlation is tighter and steeper The 
objects that present important intermediate-age population tend 
to larger AC frequency, suggesting that ACs formed from single 
stars are also present. 

The detection and characterization of ACs at much greater 
distances from the LMC centre would be very interesting. In 
fact, given the strong gradients in the LMC stellar population, 
starting from « 8 Kpc the young population completely disap- 
pears. The detection of ACs in these regions of the LMC would 
strongly point to the binary origin. 
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